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Bow shocks can be formed around planets due to their interaction with the coronal medium of the host stars. The net 
velocity of the particles impacting on the planet determines the orientation of the shock. At the Earth's orbit, the (mainly 
radial) solar wind is primarily responsible for the formation of a shock facing towards the Sun. However, for close-in 
planets that possess high Keplerian velocities and are frequently located at regions where the host star's wind is still 
accelerating, a shock may develop ahead of the planet. If the compressed material is able to absorb stellar radiation, then 
the sig nature of bow shocks may be observed during t ransits. Bow-shock models have been investigated in a series of 
papers iLlama et alll201 lllVidotto et al.Ll201(j| . l201 la lbh for known transiting systems. Once the signature of a bow-shock 
is observed, one can infer the magnetic field intensity of the transiting planet. Here, we investigate the potential to use 
this model to detect magnetic fields of (hypothetical) planets orbiting inside the habitable zone of M-dwarf stars. For 
these cases, we show, by means of radiative transfer simulations, that the detection of bow-shocks of planets surrounding 
M-dwarf stars may be more difficult than for the case of close-in giant planets orbiting solar-type stars. 
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1 Introduction 

The interaction of a planet with the corona of its host star 
can give rise to the formation of shocks that surround the 
planet's magnetosphere. Similar to what occurs around the 
Earth and other planets in the solar system, bow-shocks may 
develop around exoplanets. This idea had recently been ap- 
plied to explain the lightcurve asymmetry observed in the 
near-UV transit of the close-in giant planet WASP- 12b. 

Fossati et al. I d2010bl) observed that the near-UV transit 



lightcurve of WASP- 12b shows an early ingress when com- 
pared to its optical transit. Such observations indicate the 
presence of asymmetries in th e exosphere of the planet. In 
particular, IVidotto et al. I d2010l) suggested that this asymme- 
try could be explained by the presence of a shock surround- 
ing the planet's magnetosphere. 

To test this idea, Llama et al. ( 201 lb performed Monte 
Carlo radiation transfer simulations of the near-UV tran- 
sit of WA SP- 12b. Their results support the hypothesis pro- 
posed by Vidotto et al. ( 2010h as it explains both the ob- 
served level of absorption and the time of the (early) ingress 
observed in the near-UV lightcurve of the planet. 

An interesting outcome of this model is that one can 
probe the presence of planetary magnetic fields. Magnetic 
fields may provide an import shield against erosion caused 
by the impacting particles from the stellar wind. In the 
Earth, the presence of a magnetic field is associated with the 
development of life. M-dwarf stars are currently the main 
targets in searches for terrestrial habitable planets. There- 
fore, it is interesting to analyse whether shocks could be 
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used to trace magnetic field of Earth-like planets orbiting 
M-dwarf stars. 

The present paper is organised as follows. Section |2] 
summarizes the general characteristic of the bow-shock 
model that was initially developed for the case of WASP- 
12b. We also provide an overview of the results obtained 
through radiative transfer simulations, and discuss how this 
model could be useful to probe for planetary magnetic fields 
in other extra-solar planets. In Section[3] we investigate the 
potential for detecting bow shocks surrounding the magne- 
tospheres of hypothetical transiting planets orbiting inside 
the habitable zones of M-dwarf stars. Section|4]presents our 
conclusions. 



2 The Shock Model 

A bow shock around a planet is formed when the relative 
motion between the planet and the stellar corona/wind is 
supersonic. The shock configuration depends on the direc- 
tion of the flux of particles that arrives at the planet. We 
illustrate two different limits of the shock configuration in 
Fig.Q] where 9 is the deflection angle between the azimuthal 
direction of the planetary motion and n is a vector that de- 
fines the outward direction of the shock. As seen from the 
planet, n is the velocity of the impacting material. 

The first shock limit, a "dayside-shock", occurs when 
the dominant flux of particles impacting on the planet arises 
from the (mainly radial) wind of its host star. For instance, 
the impact of the supersonic solar wind forms a bow shock 
at the dayside of Earth's magnetosphere (i.e., at the side that 
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(a) 

d ays i de-shock 
(6=90°) 



optical wavelengths, we applied the model presented in 



Sectio nlHto the particular case of WASP- 12b dVidotto et al 
l2010h . 

WASP- 12b orbits a late-F main-sequence star with mass 
ill* = 1.35 Mq and ra dius R* = 1.57 R ^, at an orbital ra- 

120091) . Due to its close 



§4 Au~0 dius of a = 3.15 R* dHebb et al 



"r 



(b) 

ahead-shock 
(9=0°) 



Aulln 



u r ~0 



Fig. 1 Sketch of shock limits: (a) dayside-shock (8 = 
90°), (b) ahead-shock (8 = 0°). Arrows radially leaving 
the star depict the stellar wind, dashed semi-circles repre- 
sent the orbital path, 8 is the deflection angle between the 
shock normal n and the relative azimuthal ve locity of the 
planet Au. Adapted from lVidotto et alJ d201dh . 



faces the Sun). This condition is illustrated in Fig.[T^ and is 
met when 

u r > c s , (1) 

where u r and c s are the local radial stellar wind velocity and 
sound speed, respectively. 

A second shock limit, an "ahead-shock" (Fig. [TJ)), oc- 
curs when the dominant flux of particles impacting on the 
planet arises from the relative azimuthal velocity between 
the planetary orbital motion and the ambient plasma. This 
condition is especially important when the planet orbits at 
a close distance to the star, and therefore, possesses a high 
Keplerian velocity uk ■ In this case, the velocity of the par- 
ticles that the planet "sees" is supersonic if 

Au = \u K —u v \> c s , (2) 

where u v is the azimuthal velocity of the stellar corona. 

As we show later on Section 3, planets orbiting in the 
habitable zones of M-dwarf stars fall into an intermediate 
case, where both the wind and the azimuthal relative veloc- 
ities will contribute to the formation of a shock around the 
planet. 

2.1 Near-UV Early Ingress of WASP-12b 



Motivated 



dFossati et al 



b y obse rvations 
' 2010bh 



made with the HST 
which showed that in the near- 
UV, the transit of WASP- 12b starts earlier than in the 



proximity to the star, the flux of coronal particles impact- 
ing on the planet should come mainly from the azimuthal 
direction, as the planet moves at a Keplerian orbital veloc- 
ity of u K = (GA'h/a) 1 / 2 ~ 230 km s" 1 around the star. 
Therefore, stellar coronal material is compressed ahead of 
the planetary orbita l motion, possibly fo rming a bow shock 
ahead of the planet. Vidotto et al. ( 2010h , s suggestion is that 
this material is able to absorb enough stellar radiation, caus- 
ing the early-ingress observed in the near-UV light curve 
(see Figure|2]l. 

By measuring the phases at which the near-UV and the 
optical transits begin, one can derived the stand-off distance 
from the shock to the centre of the planet. In the geomet- 
rical consideration made next, we assume that the planet is 
fully superimposed on the disk of the central star, which is 
a good approximation for, e.g., small planets and transits 
with small impact parameters. Consider the sketches pre- 
sented in Fig. I2 where d op and d\jy are > respectively, the 
sky-projected distances that the planet (optical) and the sys- 
tem planet+magnetosphere (near-UV) travel from the be- 
ginning of the transit until the middle of the optical transit 

62)1/2 



Rv 



and 



(Rl-b^+TM, 



(3) 



(4) 



where b is the impact parameter derived from transit obser- 
vations, R p is the planetary radius, and tm is the distance 
from the shock nose to the center of the planet. The start of 
the optical transit occurs at phase (f>\ (point 1 in Fig. |2^), 
while the near-UV transit starts at phase <f>y (point 1' in 
Fig- Wp)- Taking the mid-transit phase at <f> = <fi m = 1, 
we note that d op is proportional to (1 — <j>i), while dyv is 
proportional to (1 — <j>i>). Using Equations (O and (01, the 
stand-off distance tm is derived from observed quantities 

1 - 



1 - 



[{Rl-b^ + R p ]-(Rl-b 



2U/2 



(5) 



We assume the stand-off distance to trace the extent of 
the planetary magnetosphere. At the magnetopause, pres- 
sure balance between the coronal total pressure and the 
planet total pressure requires that 

[B p (tm)? 



p c Au 



[B c (a) 



Pc = 



Pp, 



(6) 



8tt 1 8tt 

where p c , p c and B c (a) are the local coronal mass den- 
sity, thermal pressure, and magnetic field intensity, and p p 
and B p (tm) are the planet thermal pressure and magnetic 
field intensity at rj«. In the case of a magnetised planet, the 
planet total pressure is usually dominated by the contribu- 
tion caused by the planetary magnetic field (i.e., p p ~ 0). 
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a) Optical 
transit 




b) Near-UV 
transit 



Fig. 2 Sketches of the light curves obtained through 
observations in the a) optical and b) near-UV, where 
the bow shock surrounding the planet's magnetosphere is 
also able to a bsorb s tellar radiation. Figure adapted from 
Vidotto etal.ld2011bh . 



Vidotto et all d2010J) showed that, because WASP- 12b 
is at a close distance to the star, the kinetic term of the 
coronal plasma may be neglected in Eq. ©. They also 
neglected the thermal pressure, so that Eq. (O reduces to 
B c (a) ~ B p {tm)- Further assuming that stellar and plane- 
tary magnetic fields are dipolar, we have 



Bp — Bs, 



R* /a 
R P /r M 



(7) 



where B« and B p are the magnetic field intensities at the 
stellar and planetary surfaces, respectively. 

Therefore, by determining the phase at which the near- 
UV transit begins, one can derive the stand-off distance 
(Eq.|5]) and then estimate the intensity of the magnetic field 
of the planet (Eq. [7J, provided that the stellar magnetic field 
is known. 

For WASP- 12, we use the upper limit of <C 10 G 
dFossati et al. . l2010al) and the stand-off distance obtained 
from the near-U V transit observation tm = 4.2 R p 
( Lai et al. , 20101) . and we predict a planetary magnetic field 



2.2 Radiation Transfer Simulations of the near-UV 
transit 

In order to test the hypothesis of the bow-shock model 
was indeed able to cause the lightcurv e asymmetry ob- 
served in WASP- 1 2b. iLlama et al.1 d20 1 lb performed Monte 
Carlo radiation transfer simulations of the near-UV tran- 
sit of WASP- 12b. The characteristics of the stellar coro- 
nal plasma (d e nsity, velocity and temperature), modelled by 



Vidotto et al. ( 2010l) . were adopted in order to derive the 



density at the sh ock nose and the ang le at which the shock 



is formed. As in iVidotto et al.1 (120101) . ILlama et al.1 (1201 II) 



ngh 
10), 



assumed a shock in the adiabatic limit where the density be- 
hind such a shock increases by a factor of four with respect 
to the density ahead of the shock (stellar coronal material). 

The characteristics of the local plasma surrounding the 
planet wer e derived based on simple models that we de- 
scribe next. Vid otto et al. I (l2010h adopted two scenarios. The 
first one considers the corona as a hydrostatic medium, so 
that it corotates with the star. The other one, considers that 
the corona is filled with an expanding, isothermal wind. In 
this case, the wind radial velocity u r is derived from the 
integration of the differential equation along the radial co- 
ordinate r 

du r 1 dp c GM, 

r dr 



(8) 



p c dr r 2 

where G is the gravitational constant. The first prescription 
adopted has the advantage of having analytical solutions. 
The latter one, although lacking an analytical expression, 
can be easily integrated. In the first scenario, there is no 
radial velocity of the wind plasma, so that the shock forms 
ahead of the planet, while in the second scenario, the shock 
forms at an inter mediate angle . Thes e angles were used in 
the simulations of lLlama et al.l (1201 II) . 

To compute the near-UV lightcurve, the characterisation 
of the three-dimensional geometry of the shock is required. 
Therefore, two unknowns of the shock geometry, its solid 
angle and its thickness, had to be specifi ed. In order to tackle 
the influence of these two parameters, ILlama et al.l (1201 lb 
performed simulations for several shock geometries. As a 
result, they found that different sets of parameters could pro- 
duce similar solutions. 



To constrain the model parameter, ILlama et al.l (120111) 
rel ied on the informatio n present in the near-UV lightcurve 
bv lFossati et al.l (120 10bb . By analysing different models that 
were equally able to provide a good fit to the HST data, they 
were able to place constraints on 1) the phase of the near- 
UV ingress (<fii>) and 2) the optical depth of the shocked 
material, related to the extent of the shock and its thickness. 
From 1), they could constrain the projected stand-off dis- 
tance as 5.5 i?„, slightly la rger than the value d erived by 



Lai et al.l (1201 



-±p, angling mail my mm, u uivm u_) 

(J and used in I Vidotto etal] d2010l) . From 2) 



they showed that the shocked material does not need to have 
a large optical depth to cause the amount of absorption ob- 
served in the near-UV HST lightcurve. 



of B, n < 24 G for WASP- 12b. 



The simulations presented by ILlama et al.l (1201 lb sup- 
port the hypothesis that a bow-shock could generate an early 
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ingress of the transit, as the addition of a bow shock breaks 
the symmetry of the transit lightcurve. Nevertheless, the 
current data is not yet adequate to fully test this prediction. 
Near-UV observations of WASP- 12b (and other exoplanets, 
see next are desirable in order to test and constrain the mod- 
els. 

2.3 Magnetic Fields in Other Exoplanets 

To extend the previous model t o other transit i ng sys tems, 
near-UV data must be acquired. IVidotto et al.l (1201 lal) pre- 
sented a classification of the known transiting systems ac- 
cording to their potential for producing shocks that could 
cause observable light curve asymmetries. The data used are 
from the compilation in http://exoplanet.eu (Sept/2010) and 
the sky-projected s tellar rotation velocities u ro t sin(i) from 
Schlaufma assuming that sin(i) ~ 1. 



Once the conditions for shock formation are met, to be 
detected, it must compress the local plasma to a density high 
enough to cause an observable level of optical depth. For a 
hydrostatic, isothermal corona, the local density at the or- 
bital radius a is 

,2 



"0 



exp 



'K 



1 



a 



v 2 

n 2 



a 



1 



(9) 



where no is the number density at the base of the corona. 
The maximum coronal temperature that could still allow 
shock formation is such that c s > \u v — uk\, i.e., 

T^ax = (u v - U K ) 2 m/k B , (10) 

where the isothermal sound speed is c s = fc^T/m, fcs is 
the Boltzmann constant, and m the particle mass. If we as- 
sume that the corona corotates with the star, at the position 
of the orbital radius of a close-in planet, u v = v m ta/R*. 
By adopting the maximum temperature that could still al- 
low shock formation, we can obtain a minimum density 
required for shock formation using Equation (|9). Fig. 
shows this critical density as a function of orbital radius 
for a range of planets. This result was produced assuming 
that all stars have a base coronal density equal to that of 
the Sun (10 8 cm -3 ), which is a reasonable approximation 
for stars with spectral types similar to that of our Sun. The 
most promising candidates to present shocks are: WASP- 
19b, WASP-4b, WASP- 18b, CoRoT-7b, HAT-P-7b, CoRoT- 
lb, TrES-3, and WASP-5b. 



Vidotto et al.l d2011al) estimated that a lower limit of the 
stellar coronal density n m ; n ~ 10 4 cm~ 3 could still pro- 
vide detection (dashed line). Taking that into account, we 
note that a reasonable number of planets should lie above 
this density-detection threshold, suggesting that a detectable 
shock might be a common feature surrounding transiting 
planets. 

3 Shock Formation around Planets Orbiting 
M-dwarf Stars 

Current and future missions like Kepler, CoRoT and Plato 
are designed to search for Earth-like planets. Low-mass 




0.02 0.03 0.04 0.05 0.06 0.07 

Orbital radius [AU] 

Fig. 3 Predicted coronal density at the orbit of each tran- 
sit planet known as of Sept/2010, assuming that the stellar 
magnetic field is strong enough to confine the coronal gas 
out to the orbit of the planet. The dashed line represents 
a lo wer limit for detection of bow shocks. Figure adapted 
from lVidotto et all d2011al) 



stars (~ 0.06-0.8M Q ) are currently the main targets in 
searches for terrestrial habitable planets as their low bright- 
ness provide good contrast for detection of smaller-radius 
planets (transit method) and their low masses facilitate 
detection of lower-mass planets (radial velocity method). 
These methods favour detection of planets orbiting close to 
their host star, which in the case of low-mass stars is also 
the region where the habitable zone (HZ) is located. The 
concept of habitable zones relies on the location where liq- 
uid water may be found, which is essentially determined 
by the luminosity of the star. Because M-dwarf stars are 
low-luminosity objects, the habitable zone is believed to lie 
near the star. The habitable zone for a 0.5M Q star is be- 
lieved to lie between 0.2 - 0.5 AU, while for a 0.1M Q star, 
it is thought to lie m uch closer, between 0.02 - 0.05 AU 



( iKasting et al.Ul993l) . 



However, even if a planet is in the HZ of a low-mass 
star, the interaction of the planet with the host star's wind 
may erode its atmosphere, removing an important shield for 
the creation and development of life. The interaction of a 
planet with its host star wind can produce a shock surround- 
ing the planet. As presented before, if we are able to detect 
the stand-off distance between the shock and a planet (e.g., 
as during near-UV transit observations), we are then able to 
estimate the planet's magnetic field intensity. 

In this section, we calculate the characteristics of the 
wind of two M-dwarf stars with spectral types Ml. 5 and 
M4.0, adopti ng a thermally driven wind (Equation ((H), 



ParkerL 119581) . We assume these winds are isothermal with 
temperature T = 1 X 10 6 K and their densities are such 
that their mass-loss rates are M = 10 -13 M Q yr~ J . TableQ] 
shows the characteristics adopted for the star and the planet 
in our analysis. The top panels in Figure |4] show the wind 
velocity solution (solid line) and the density profile (dashed 
line) for the Ml. 5 case (left) and the M4.0 case (right). 
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Fig. 4 Results obtained from our wind models for the Ml. 5 case (left) and M4.0 case (right). Top panels: wind velocity 
solution (solid line) and the density profile (dashed line). Middle panels: magnetospheric size as a function of orbital radius 
assuming a planet with a magnetic field similar to Jupiter (B p = 14 G, dashed line) and with a magnetic field similar to 
Earth (B p = 1 G, solid line). Bottom panels: the shock angle as a function of the orbital radius of a planet. 
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Table 1 Values adopted in our models (upper part) and 
results obtained (lower part). 



Quantity 




Ml. 5 


M4.0 


M*(M ) 




0.57 


0.28 


R*(R@) 




0.51 


0.34 


a(AU) 




0.2 


0.03 


en 




82 


71 


(r M /Rp)[B p = 


= 14 G] 


12.2 


6.7 


(r M /Rp)[B p = 


= 1 G] 


5.1 


2.8 


"Mg, shock (cm 


- 3 ) 


0.09 


5.6 



a = 0.2 AU orbiting the Ml. 5 case, this angle is 6 = 82°, 
while for a planet at a = 0.2 AU orbiting the M4.0 case, 
this angle is 8 = 71°. We note that, due to the geometry of 
the transiting systems, dayside-shocks are more difficult to 
be detected as, during transit, the planet occults the shock 
from the observer. Furthermore, this type of shock may not 
be able to brake the symmetry of the transit lightcurve. Nev- 
ertheless, for the hypothetical planet orbiting inside the HZ 
of the M4.0 star, the angle we found (6 = 71°) could still 
produce an early ingress. 



3.1 Deriving the Magnetospheric Characteristics 

In the calculations presented next, we assume the planet to 
be in a circular orbit. We use Equation (|6]l to derive the plan- 
etocentric distance to the magnetopause, such that 



3.2 Detecting Bow-Shocks 



_3_ 

87tR, 



1/6 



(ID 



where we prescribed the planetary magnetic field topology 
to be dipolar. The total pressure exerted in the planet due to 
the coronal wind P c ft is 

\B c {ar 



P off = p c Au 



Pc 



(12) 



Our wind solutions neglect the presence of the stellar mag- 
netic field {B c = 0). We note that the local (i.e., at r = a) 
coronal wind velocity u r and density p c are obtained from 
the wind models. In the case of WASP-12b, data allowed 
us to estimate Tm and then derive B p . Here, we take the 
inverse procedure. 

The middle panels in Figure 0] shows the magneto- 
spheric size the planets orbiting the M-dwarf stars anal- 
ysed here would present as a function of orbital radius. We 
considered two different cases and assumed a planet with a 
magnetic field similar to Jupiter (B p = 14 G, dashed lines) 
and with a magnetic field similar to Earth (B p = 1 G, solid 
lines). For the Ml. 5 case, a planet orbiting at a = 0.2 AU 
(in its HZ) would present a magnetospheric size of tm = 
12.2i? p and 5.1 R p for B p = 1 G and 14 G, respectively. 
For the M4.0 case, the same hypothetical planet orbiting at 
a = 0.03 AU (in its HZ) would present tm = 6.7 R p or 
2.8 R p for B p = 1 G and 14 G, respectively. 

The orientation of the shock normal depends not only 
on the wind velocity, but also on the Keplerian velocity of 
the planet and on the azimuthal velocity of the wind. In the 
Parker wind solution, rotation is neglected. So the cases we 
study here are valid for slowly rotating stars. This angle is 
given by 



/ u r 
atari 



(13) 



and is 8 = 90° for a dayside-shock and 8 = 0° for an ahead- 
shock (Figure Q]). The bottom panels in figure [4] shows the 
shock angle as a function of the orbital radius of a planet. 
As we can see, the farther out the planet is, the shock 
tends to approach the dayside-type shock. For a planet at 



As discussed in Section l23l the density of the shocked ma- 
terial should be high enough as to produce a significant level 
of absorption for it to be detected. Using the density derived 
in our wind models, we calculate the density in magne- 
sium assuming the stars have solar metallicityU Because the 
shock compresses material ahead of the planetary motion, 
the density of the shocked material increases by a certain 
factor. For a strong shock in the adiabatic limit, this increase 
has its maximum value of 4. Using the results for the wind 
models we derived, we estimate this increase to be of a fac- 
tor of 3 for both stars analysed here. Therefore, the densities 
of the shocked material for both cases are shown in TableQ] 
As it can be seen, these values are considerably smaller 
than the value deriv ed for WASP-12b n Mg ~ 400 cirT 3 
d Vidotto et all 1201 Oh . indicating that the detection of bow- 
shocks of planets surrounding M-dwarf stars may be more 
difficult. 

Using radiativ e transfer s i mulati ons similar to the 
ones presented in iLlama et alj (1201 lb . we computed the 



lightcurve of a planet orbiting in the HZ of the stars mod- 
eled here (values are shown in Table [TJ. With reasonable 
assumptions for the shock extension and thickness, a bow- 
shock surrounding a 2 i?Earth -planet can not be detected, as 
it does not present a sufficient amount of optical depth that 
causes a detectable absorption. 

It is worth noting that our results depend on the con- 
ditions adopted for the wind. For instance, mass-loss rates 
in M-dwarf stars are poorly constrained. A higher value of 
mass-loss rate would result in a higher value of the wind 
density, and therefore increased density in the shocked ma- 
terial that could be detected during transit. To illustrate that, 
we computed the lightcurve of a 2 i?Ncptunc-planet orbiting 
a M4.0-star at an orbital radius a = 0.03 AU. We adopted 
a density for the wind (and for the shock) that is one or- 
der of magnitude larger than the one shown in Table [T] im- 
plying that such a star would present a mass-loss rate of 
M = 10~ 12 M Q yr _1 . We assumed that the planetocentric 
distance to the magnetopause is tm — 5 R p (or equiva- 
lently, B p ~ 18 G), the shock thickness is 1.5 R p and its 
area lies within a solid angle of A6 = 50°. The shock nor- 
mal is oriented at 8 = 71°. Figure [5] shows a comparison 



We calculate the density of magnesium as a proxy for the detection of 
the shock in the near-UV, in analogy to the near-UV transit observed for 



WASP- 12b dFossati et alll2010bl) . 
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of two synthetic lightcurves: one case where no bow shock 
is present (dashed line) and one case with bow shock (solid 
line). Note that our choice of parameters for the fictitious 
planet leads to a transit depth of ~ 5 percent normalised 
to the out-of-transit flux for the case without a bow shock 
and a ~ 7 percent-depth for the case where the bow shock 
was considered in the simulations. We note that the early- 
ingress in t his case is less-pron ounced than in the case of 



WASP-12b (Llama et aU 120111) . 




Without bow shock 
With bow shock 



wind. We showed that shock can be formed surrounding the 
planetary magnetosphere. For the case of a planet orbiting 
at 0.2 AU from the Ml. 5 star, we showed that the shock 
forms at an angle closer to 9 = 90°, implying that the shock 
should be occulted from the observer during transit. If this 
is the case, the presence of the shock may not brake the op- 
tical lightcurve symmetry. For the case of a planet orbiting 
at 0.03 AU from the M4.0 star, the shock forms at a smaller 
angle = 71°. 

However, for the shock to be detected by means of an 
early-ingress, it has to compress the stellar coronal material 
to a level of optical depth that can cause enough absorption 
during transit. For the wind characteristics we assumed, the 
density we computed in the shocked material does not seem 
to be significantly high to cause an observable early-ingress 
in the planet near-UV lightcurve. This result was confirmed 
by means of radiative transfer simulations. This implies that 
the method applied to close-in giant planets orbiting solar- 
type stars may not be as useful in detecting magnetic fields 
of planets orbiting inside the habitable zones of M-dwarf 
stars. 



Fig. 5 Lightcurves of a 2 i?Neptunc -planet orbiting a 
M4.0-star at an orbital radius a = 0.03 AU. Dashed line 
shows the lightcurve produced when no bow-shock sur- 
rounds the planet and the solid line is the result when the 
bow shock is included in the Monte Carlo radiation transfer 
simulation. 



We also note that M dwarf stars can exhibit substantial 
activity, which can lead to variations in the ambient medium 
surrounding the planet. In response to these variations, the 
planet's magnetosphere will adjust itself. Because of that, 
transit lightcurves in the near-UV can present variability 



(IVidottoet all 1201 lbh 



4 CONCLUSIONS 

In this paper, we investigated the potential for detecting bow 
shocks surrounding the magnetospheres of transiting plan- 
ets orbiting inside the habitable zones of M-dwarf stars. 
The idea that bow shocks may be present arou nd exoplanets 



have now been developed i n a series of papers (ILlama et al 



201 li IVidotto et all l2010ll201 lalbl) . Due to the interaction 
of the planet with the wind of its host star, bow shocks may 
form surrounding the magnetosphere of planets, and, if ob- 
servable, the shock can help constrain planetary magnetic 
fields. 

In order to analyse whether shock could be formed on 
planet s orbiting M-dwarf stars, we applied IVidotto et al 
(2010)'s model to two hypothetical low-mass stars with 



spectral types Ml. 5 and M4.0. Earth-size planets were as- 
sumed to orbit inside the habitable zones of their host-stars 
and the stars were assumed to have a thermally driven-type 
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